The Back and Forth Nudging algorithm for data 
assimilation problems: theoretical results on 
transport equations 

Didier Auroux*^ Maelle Nodet ^§ 
September 29, 2008 

Abstract 

English: 

In this paper, we consider the back and forth nudging algorithm that has 
been introduced for data assimilation purposes. It consists of iteratively 
and alternately solving forward and backward in time the model equa- 
tion, with a feedback term to the observations. We consider the case 
of 1-dimensional transport equations, either viscous or inviscid, linear or 
not (Biirgers' equation). Our aim is to prove some theoretical results on 
the convergence, and convergence properties, of this algorithm. We show 
that for non viscous equations (both linear transport and Biirgers), the 
convergence of the algorithm holds under observability conditions. Con- 
vergence can also be proven for viscous linear transport equations under 
some strong hypothesis, but not for viscous Biirgers' equation. Moreover, 
the convergence rate is always exponential in time. We also notice that the 
forward and backward system of equations is well posed when no nudging 
term is considered. 

French: 

Ce travail etudie I'algorithme du nudging direct et retrograde, qui a ete 
introduit en assimilation de donnees. Get algorithme consiste a resoudre 
iterativement I'equation du modele, agrementee d'un terme de rappel aux 
observations, dans le sens direct puis dans le sens retrograde. Dans ce 
travail nous nous interessons aux equations de transport en dimension 
1, avec ou sans viscosite, lineaires ou non (Biirgers). Notre objectif 
est d'etudier la convergence eventuelle, et la vitesse de convergence le 
cas echeant, de cet algorithme. Nous prouvons que, pour les equations 
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non visqueuses (lineaire ou Burgers), la convergence a lieu sous des hy- 
potheses d'observabilite. La convergence peut aussi etre demontree pour 
des equations de transport lineaires visqueuses sous des hypotheses fortes, 
mais pas pour I'equation de Biirgers visqueuse. En outre, lorsque la 
convergence a lieu, la vitesse de convergence est toujours exponentielle 
en temps. Nous remarquons aussi que le systeme d'equations directe 
et retrograde est toujours bien pose lorsqu'aucun terme de rappel n'est 
present. 

1 Introduction and main results 

Data assimilation is the set of techniques aiming to combine in an optimal way 
the mathematical information provided by the model equations and the physi- 
cal information given by observations, in order to retrieve the state of a system. 
Several types of methods have been widely studied in the past decades. We 
can cite here interpolation, variational and stochastic methods. The first ones 
interpolate the measurements from the points of observation towards the grid 
points, the interpolation being weighted by the statistics of the observations 
|Kal03j . Variational methods are based on the optimal control theory, and data 
assimilation is set as being a problem of constrained optimization. The goal is 
to minimize a cost function measuring the difference between the observations 
and the corresponding quantities provided by a model integration. The initial 
condition of the system can then be seen as a control vector [LDT86j . Finally, 
the basic idea of stochastic methods is to consider the fields as the realization 
of a stochastic process and carry out Kalman filtering methods [Kal60[ IEvL99| . 
We can also mention one of the very first data assimilation schemes: the nudg- 
ing method. Also known as Newtonian relaxation or dynamic initialization, it 
consists of adding a feedback term to the observations directly in the model 
equations [HA76| . 

All these methods require extensive work, either from the implementation 
or the computation point of view. For instance, variational methods require 
the linearization of all operators and also the implementation of the adjoint 
model. They also need efficient optimization schemes, as the minimization is 
performed on spaces of huge dimension. On the other side, stochastic methods 
are somewhat easier to implement, but they require knowledge, storage and 
manipulations of huge matrices. 

The Back and Forth Nudging (BFN) algorithm has recently been introduced 
as a simple and efficient method for solving data assimilation problems [AB05] . 
In most geophysical applications, data assimilation consists of estimating a tra- 
jectory, solution of a partial differential equation (PDF), from the knowledge 
of observations. These observations are usually sparse in time and space, and 
incorrect in the sense that they are not the restriction of a solution of the PDF 
model. One step of the BFN algorithm consists of solving first the model equa- 
tion, in which a feedback to the observation solution is added, and then the same 
equation but backwards in time, with also a feedback term to the observations. 
Such forward and backward integrations provide a new value of the solution at 
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the initial time t = and the aim of the BFN is to improve the quahty of the 
initial condition. 

The idea of the back and forth nudging is to use the difference between 
the observations and the model trajectory as a feedback control of the equa- 
tions, both in the forward and backward integrations. This makes the numerical 
scheme extremely easy to implement, in comparison with both variational and 
stochastic methods, as we usually only consider diagonal (or even scalar) gain 
matrices. The back and forth nudging scheme can also be seen as an intermedi- 
ate scheme between variational and stochastic methods, as the standard nudging 
technique has both variational (minimization of a compromise between the ob- 
servations and the energy of the system) and stochastic (sub-optimal Kalman 
filter) interpretations |AB08| . 

As a first approximation, we consider in this paper that the observations 
are correct (i.e. no observation error), and hence the observations satisfy the 
model equation. We consider various observation domains: first we assume that 
the observations Uobs{t,x) are available for any point x and time t, second we 
assume that they are available for t S [ti, t2] and for all x, and third we consider 
that they are available for all t over a given space domain. This is done through 
the time and space dependency of the feedback (or nudging) gain matrix K(t, x) 
that is equal to when the observations are not available. 

Many numerical experiments in almost realistic situations suggest that this 
algorithm works well, and that the identified solution gets closer to the ob- 
servations [AB08| . The goal of this paper is to prove some theoretical results 
and convergence properties in the particular case of transport equations, either 
viscous or inviscid, either linear or non-linear (Biirgers' equation). 

In section [21 we consider one step of the BFN algorithm applied to a linear 
viscous transport equation: 



I u\t=o = "0 

(1) 

{dtu - vdxxU + a{x)dxU K'{u-Uobs) 
u|x=o = u\x=i = 
u\t=T = u{T) 

where the following notations hold for all further cases: 

• the time period considered here is i e [0, T]; 

• the first equation [F) is called the forward equation, the second one {B) 
is called the backward one; 

• K and K' are positive and may depend on t and x, but for simplicity 
reasons, we will always assume that there exists a constant k G M^j. such 
that K'(t,x) = KK{t,x); 

• a{x) € W^''°°{fl), ri being the considered space domain, either the interval 
[0, 1] or the torus [0, 1]; 








K{u 
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• 1/ is a constant; 

• Uobs is a solution of the forward equation with initial condition Uo^^: 



dtUobs - i^d^xUobs + a{x)dxUobs = 

mU=o = u\x=i (2) 

''obs 



u\t=o = u° 



The following result holds true: 



Theorem 1 We consider the one-step BFN (QP with observations Uobs satisfy- 
ing We denote 

W{t) = u{t)-Uobs{t) /o'l 
W{t) = u{t)~Uobs{t) ^ ' 

Then we have: 

1. If Kit, x) ~ K, then we have, for all t G [0,T]; 

w{t) = e(-^'-^')(^-*)iz;(i) (4) 

2. If K(t,x) = K{x), with Support (K) C [a,b] where a < b and a ^ or 
b 1, then equation (Qp is ill-posed: there does not exist a solution {u,u), 
in general. 

3. If K{t,x) = Klyt^^t^^it) with <ti <t2 <T, then we have 

^5(0) = e(-"')(*^-*i)w;(0) (5) 

In section [U we consider one step of the BFN algorithm applied to the 
viscous Biirgers' equation: 

dfU - vdxxU + ud^u = -K{u - Uobs) 
(F) { mU=o = '«U=i = 

(6) 

dtu - vdxxU + udxU = R'{u-Uobs) 
(B) { u\x=o = u\x=i = 

u\t=T = u{T) 

with the same notations as before. 

The observations Uobs satisfy the forward Burgers' equation: 

dtUobs - vdxxUobs + UobsdxUobs = 

u\x=o=u\x=i = (7) 



We have the following result ii K ^ 0: 
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Theorem 2 The BFN iteration for viscous Biirgers' equation, with obser- 
vations Uobs satisfying is ill-posed, even when K{t,x) is constant (except 
for K{t,x) = Q): there does not exist, in general, a solution {u,u). 

In the particular case when A' = A'' = 0, the backward problem is ill-posed 
in the sense of Hadamard, but it has a unique solution if the final condition u\t=T 
is set to a final solution of the direct equation. Moreover, in this particular case, 
the backward solution is exactly equal to the forward one: u{t) = u[t) for all 
t G [0,T]. The main result is the following: 

Proposition 1 If K ^ K' = 0, then problem is well-posed in the sense of 
Hadamard, and there exists a unique solution {u,u). Moreover u~u. 

Section [4] considers the extension of theorem [1] to the inviscid case, for both 
linear transport and Biirgers' equations. 
We first consider the linear case. The BFN equations are: 



{F) { 



'JtU 



a{x)dxU 
u\x=o 
dxu\x=o 

a{x)dxU 
u\x=o 

dxu\x=o 
u\t=T 



-K{u- Uobs) 

u\x=l 

dxU\x=l 

Uq 

K'{u - Uobs) 

u\x=l 

dxU\x=l 

u{T) 



(8) 



where a{x) can be constant or not. 

Theorem 3 We consider the non viscous one-step BFN with observations 
Uobs satisfying with K = Q. We denote 



wit) 
wit) 



uit) - Uobsit) 
Uit) - Uobsit) 



(9) 



We denote by 

is,^is,x)) (10) 

the characteristic curve of equation ^F) with K = 0, with foot x in time s — Q, 
i.e. such that 

(s,7/.(s,x))|,,=o = (0,x) (11) 

We assume that the final time T is such that the characteristics are well defined 
and do not intersect over [0,r]. 
Then we have: 



1. If Kit, x) = K, then we have, for all t G [0,T]; 

wit) = u;(t)e(-^-^')(^-*) 



(12) 
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2. If K{t,x) = Kl^tiM]i't) '^^ih < ti < t2 < T , then we have 

w(0) = w(0)e(-^-^')(*^-*i) (13) 

3. If K(t,x) — K{x), then we have, for all t G [0,T]; 

w{t, jpit, x)) = w{t, -tpit, x)) exp ^- K{->p{s, x)) + K'{ip{s, x)) ds^ 

(14) 

We finally consider non viscous Biirgers' equation, still with periodic bound- 
ary conditions, and for a time T such that there is no shock in the interval 
[0,T]: 

dfU + udxU = -K{u-Uobs) 



u\x=0 = u\x=l 

dxu\x=o = dxu\^^i 

u\t=o ^ Uq 

(15) 

dtu + udxU = R'{u-Uobs) 

u\x=0 = u\x=l 

dxu\x=o = dxu\x=i 

u\t=T = u{T) 



Theorem 4 We consider the non viscous one-step BFN \15]) . with observations 
Uobs satisfying \15V F) with K = Q. We denote 

w{t) = u{t)-Uobs{t) 

wit) = u{t)-Uobs{t) ^ ' 

We assume that Uobs G M^^'°°([0, T] x Q), i.e. there exists Af > such that 

\dxUobs{t,x)\ < M, Vt e [0,T],Vx e n (17) 

Then we have: 

1. If K{t,x) = K, then we have, for all t G [0,T]; 

\\w{t)\\ < e^~"~"'^^'^^^~'^Mt)\\ (18) 

2. If K{t,x) = ^^'^ < ii < i2 < r, then we have 

||w(0)|| < e(-^-^')(*^-*i)+^^^||w(0)|| (19) 

Proposition 2 We consider one forward (resp. backward) BFN step of the non 
viscous Biirgers equation il5lF) (resp. il5lB) ). With the notations of theorem 
if K{t,x) = K{x), then we have 



w(T,ip{T,x)) ^ w{0,x)cxp \ - I K{tp{a,x))da - I dxUobs{o;ip{cr,x))da j 

(20) 
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Remark 1 For the special case K{t,x) = K{x) = K\]^g^^{x) where K is a 
constant and [a, b] is a non-empty sub-interval of [0, 1], we have 



w{T,ii){T,x)) = w(0,a;)cxp I -Kxix) - I dxUobs{cr,ij{(r,x))da (21) 



where ^ 

X{x)= I l[,,b](V'(fT,a;))da (22) 
Jo 

is the time during which the characteristic curve ip{a,x) with foot x of equation 
hlSV F) with K = Q lies in the the support of K . The system is then observable if 
and only if the function x has a non-zero lower bound, i.e. m := minxix) > 0, 

the observability being defined by (see \Rus78l ): 

3Cyu solution of IMF) with K = 0, \\u{T,.)\\^<c[ \\K{.)u{s, .)\\^ ds 

In this case, proposition proves the global exponential decrease of the error, 

MT 

provided K is larger than , where M is defined by equation |i7l). 

TO ' — ' 

From remark [l] we can easily deduce that if for each iteration, both in the 
forward and backward integrations, the observability condition is satisfied, then 
the algorithm converges. Note that this is not a necessary condition, as even if 
x{x) = 0, the last exponential of equation ((2T|) is bounded. 

Note also that in real geophysical applications (either meteorology or oceanog- 
raphy), there is usually no viscosity. In this case, assuming the observability 
condition, the BFN algorithm is well posed, and theorem 3] and proposition [2] 
say that the solution tends to the observation trajectory everywhere, and not 
only on the support of K. From a numerical point of view, we can observe 
that even with discrete and sparse observations in space, the numerical solution 
is corrected everywhere jAB08| . We also observed that with a not too large 
viscosity coefficient, the behavior of the algorithm remains unchanged. 

Figure [1] illustrates the results given in theorem [3] in the case 3 (top) and 
proposition [2] and remark [1] (bottom) . These numerical results correspond to 
a simple case: Uots = 0, uo{x) ~ asin(27ra;), K = K' = l[o;o.5](^)- Various 
final times T are considered, from 0.05 to 1, and both figures show the following 
expression 

->»!(^) P3) 



^w(0, x) 

as a function of x e [0;1]. Figure [T]-top illustrates equation p^ . The best 
possible decrease rate is then max(/ir + K') x T = 2T. In the linear case, the 
transport is a{x) = 1. As half of the domain is observed, the observability 
condition is satisfied iff' T > 0.5, and this is confirmed by the figure. Concerning 
Biirgers' equation, figure [l]-bottom illustrates equation ((20)l . After one iteration 
of BFN, the best possible decrease rate is also 2T. We can see that in this 
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Figure 1: Decrease rate of the error after one iteration of BFN (see equation [23]) 
as a function of x, for various times T; top: linear transport equation; bottom: 
inviscid Biirgers' equation. 
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case, due to the nonlinearities of the model, the solution is less corrected on 
[0; 0.1] but more on [0.5; 0.6]. From this figure, we can see that the observability 
condition is satisfied for T larger than approximately 1. 
Finally, some conclusions arc given in section [S] 

2 Linear transport equation with a viscous term 

In this section we prove theorem [TJ 



2.1 Case 1: K constant 

The differences w and w satisfy the following equations: 
dtw — vdxxW + a{x)dxW + Kw 

dfW — vdxxW + a{x)dxW ~ K'w 

W\x=li = W\x=l 
w\t=T 






Wo 




w{T) 



(24) 



Wc denote by S+ and S- the operators associated to these equations, seen as 
forward equations on [to,t] with initial conditions given in U): 



S+{to,t){wito))^w{t), S^{to,t)iwito)) ^ w{t) 
The BFN algorithm has a solution if and only if wc have 

w{T) G Im(5_(0,T)) 
We re-write equation (P5|) associated to w. 



(25) 
(26) 

(27) 
(28) 



If we assume that the expected result is true, i.e. w{T) e Im(5_(0,T)), then 
we can assume that it is also true for all t, i.e. we can assume that: 



dtW — vdxxW + a(x)dxW — K'w 
w\x=o = w\x=\ 

so that we have, thanks to Duhamel's formula: 



{^-~K-K')w 




w(t) = S'_(0,<)(wo) + / S^{s,t)[ly-K - K')w{s))ds 



yt,3ip{t),w{t) = S-{0,t)ip{t) 
In that case, we replace ([^^ in ([25)1 and we get: 



(29) 



w(t) = 5_(0,t)(wo) + / S-is,t){{-K - K')S-iO,s)ip{s))ds (30) 
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As the equation is linear, the scalar coefRcient {K ^ K') commutes with S- and 
we get: 



w{t) 
So that we have 



JO ^ 

S-{0,t){wo) + {-K - K')S-{0,t) I Lp{s)ds 



iwo) + {-K-K') / ^is)ds 



ip{t) = Wo + {-K - K') / ip{s) ds 



i.e., (fi satisfies 
and finally 

so that we get for w(r) 



u;(r) = 5_(0,r)(u;oe(-^-^')^) 



(31) 
(32) 

(33) 
(34) 



(35) 
(36) 



Reciprocally, setting 

leads to w satisfying w{T) ~ w{T), so that {w, w) is the solution of the one-step 
BFN pi]) . 

Moreover, we have, for all t £ [0, T]: 



w{t) ^ S'_(0,t)(w(0)) 

= e(-^-^')(^-*)5_(0,t)(«;oe(-^-^')*) 

= e(-^-^')(^-*)5_(0,t)(</.(t)) 

= e(-^-^')(^-*)zi;(i) 



(37) 



2.2 Case 2: K{x) 

Wc assume that Support (K) C [a, b] where a < b and a Q or b 1, i.e. the 
support of K is not [0, 1]. Wc can follow the same reasoning as previously up 
to equation pO]) : 



w{t) = S^{0,t)ip{t) = S-{0,t){wo)+ I S'-(s,t) [{~K{x) ~ K'{x))S-{0,s)ipis)] ds 

(38) 

Let us assume, by contradiction, that —K(x) — K'{x) commutes with S'_. Then 
we get: 



S^{Q,t){ip{t) -wq) ^ {~K{x) - K'{x))S^{{),t) / Lp{s)ds 



(39) 
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But we know that has the unique continuation property, that is: 



Proposition 3 If S-{0,t){X) = on a non-empty subset of [0,1], then S-{Q,t){X) 
on [0,1]. 



(40) 



This resuh and pQ]) give: 

w(t) =5_(0,t)((p(t)) =5_(0,t)(wo) =5+(0,i)(wo) 



As this stands for every wq, we have S- = 5+ and finally K = K' = 0, which 
is a contradiction. Therefore, K + K' does not commute with 5_. Thus, in 
general, we cannot find any function -0 such that: 



S- (s, t) [{-K{x) - K'{x))S- (0, s)^(s)] ds = 5_ (0, t)^' (41) 



2.3 Case 3: K{t) 

We assume that K{t,x) = /-5:(t) = Kl[t,^t2]it) with < ti < t2 < T. We can 
follow the same reasoning as for K constant, up to the Duhamcl formula (|30p : 



wit) = S.iO,t)ipit) - S.iO,t)iwo)+ / 5-(s,t) [(-i^-X')l[ti.t2](s)^-(0,s)^(s)] ds 

Jo 

(42) 

As + ii"' is independent of x, it commutes with and wc have for tp: 



( 



if t < tl 



S^{0,t)ipit) = S^iO,t){wo)+{ 



{-K - K')S-{0,t) [ ip{s)ds ifti<t<t2 



{-K - K')S-{0,t) / (/?(s)ds ift<t2 



So that the corresponding is given by: 



Wo if t < tl 

^{t) = { woe(-^-^')(*-*i) ifti<t<t2 
woe(-^-^')(*=-*i) ift<t2 



(43) 



(44) 



And thus the result follows. 



3 Burgers' equation with a viscous term 
3.1 Proof of theorem [2] 

Without loss of generality we assume that the observations are identically zero: 
Uobs{t,x) = for all (t,a;). Let us first introduce some notations. 
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Let us denote by w (resp. w) the differences between u (resp. u) and the 
observations, as in ([3]), they satisfy the foUowing equations: 



dtw — vdxxW + wdxW + A'w = 

{F) { w\x=o = w\x=i = 

w\t=o = Wo 

dtW — vdxxW + wdxW — K'w = 

(B) u5U=o = w|3;=i = 

w\t=T = w{T) 



(45) 



Let us denote also by S+ and S- the non-hnear operator associated to the 
forward equations with K or K': 

S+(tQ,t){w{to))^w(t), S-{to,t)iw{to))^w{t), yt>ta (46) 

We will also use the linear operators [/+ and t/_ associated to the following 
linear equations: 

dt4>-vdxx4' + Kcj) = /n A m+\ ia'7\ 

^U.o = 0U^i = o, = 00 ^ u^{0MM = m (47) 

= = 0, 0Uo = 00 ^ C/-(O,t)(0o) = 0(t) (48) 

To prove theorem [5] we will prove that w is not in the image of 5_ , in 
general. To do so we will use perturbations theory. We can easily show that 5*+ 
is infinitely continuous with respect to the data wq . So if we suppose that wq is 
small: 

Wo = £<<5o (49) 

then we have that w(t), solution of the forward equation (^SlF) is also small 
and can be developed in series of e 



e"u." (50) 

Similarly, we develop w in series of e 



w — e 

n>0 



(51) 



e 

n>0 



As previously, w satisfies: 



dtW — vdxxW + Kw = —wdxW 

w\x=o = w\x=i = (52) 
w|t=o = Wo 



so that if we develop in series of e we get, for w^: 



dtw° - iydxxw° + Kw° = 

w''\x=o ^ w°\x=i = (53) 
w°\t=o = Lpo 
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For we have: 



dtw^ — vdxxW^ + Kw"^ = —w^dxw'^ 

w^\x=o = w^\x=i = (54) 
w^lt^o = 

Similarly we have for and w^: 

dtw° - vdxxw'^ - K'w^ = 

w°\x=o = w"\.=i = (55) 
w°\t=T = w^{T) 

dtw^ — vdxxW^ — K'w^ = —w^dxVJ^ 

w'^\x=o = w^\x=i = (56) 

We can compute w'^ and thanks to [/+: 
w^it) = C/+(0,t)(^o) 

w^it) ^ - [ U+is,t)[w"is)dxw''is)]ds ^^^^ 
Jo 

If we assume that is well defined, with 

w"{t)^U-{Q,t){i^o) (58) 
then the condition w{T) = 'w{T) leads to 

[/_(0,T)(^o) = c/+(o,r)(^o) 
=> ^0 = C/-(0,T)"i[/+(0,T)(^o) (59) 

^ ^0 = c-(^^+^')^^o 



Then we have for w^: 



d^{t) = f/_(0,i)(V'o) 

= C/_(0,i)e-(^+^')^^o 



For u)^ the final condition w^{T) = w-^{T) gives, thanks to ([57)1 : 
wi(T) = u;i(r) 

T 

L/+(s,t)[u;°(s)9:,?i;"(s)]ds 



(60) 



(61) 



On the other hand, if we assume that is well defined, with w^{0) = ipr, then 
equation ([5S)) and the Duhamel formula give 

w^{T) = U^{0,T)I^Pt]~ [ U^{s,T)[w°{s)dxW°{s)]ds (62) 



13 



Then, equalling (|62p and (pT|) wc should have 

U.{Q,T)[^t]- I C/-(s,r)[SiO(s)9,zIi°(s)]rfs = - / U+{s,T)[w\s)d^w''{s)]ds 
Jo Jo 

(63) 

Therefore 



2C/_(0,T)[V't] = r U-{s,T)[d^,{w'>isf)]ds - r U+{s,T)[d,{w°isf)]ds 
Jo Jo 

^ (64) 

If we assume that ■i/'T = ^'^xgT: then we obtain, up to a constant 



c/-(o,r)[g3 



U_{s,T)[w"{sf]ds~ r U+{s,T)[w°{s)^]ds (65) 



We now use dST]), dSS]) and ((55)) : 



[/_(s,r)[([/_(0,s)(e-(^+^')^V'o)]'ds 
[/+(s,T)[([/_(0,s)(v5o)]2ds 



(66) 



And if Js' > and K' > this last equation is in general impossible: such gx 
does not, in general, exist. Indeed, let us do an explicit computation thanks to 
Fourier series: 

(^o = $^a„e™^ 5T = ^&„e™^ (67) 



n>l 



n>l 



We recall that we have 



U+is,t) 
U-is,t) 



E 

n>l 



n>l 



(68) 



gma:g(i^' — i/n^) (i — s) 



n>l 
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Then we can compute the right hand side of equation (|66p : 



2(M')T_i) / u+{s,mU_{0,s)icp„)]^ds 



g-2(M')T_i) r u+is,T) 

Jo 

e-2(M')T_i) r U+{s,T) 



ds 



E 



JlsK' Anx 



Q-2(K+K')T -I) I 

Jo 



ds 



e e apttqB 

p+q=n 

n p+q—n 

EY^ -KT-vTri^+inx 2sK' +sK+usn'^-sy{p^+q^) 

n p+g— n 

J2TK'+TK+2upqT -r 



ds 
ds 



E E 



iCoe 



n p+g— n 

For the left hand side of (1551) we have: 



2if' + if + 2^-^(7 



C/_(0, T)[5t] = ^6„e™^e^' 

n 

So that we get, for all n: 



T-un^T 



(69) 
(70) 



6„ = e(-^'+^)^(e-2(^+^')^ - 1) 



E 



.P+9=" 



2/\ ' + if + 2vpq 



(71) 

This defines a distribution iff 6„ has polynomial growth, iff 6„ has polynomial 
growth, where 



6„ = (e-2(^+^')^ - 1) 



q2TK' +TK+2ypqT 

E 2X' + K + 2vpq 

,p+q=n 



(72) 



which is clearly not the case for every sequence (a„) with polynomial growth, 
unless K = K' ^0. 

3.2 Particular case: K = K' = 

We consider the particular case where K = K' = 0, i.e. there is no nudging 
term either in the forward or backward equations. In this case, proposition [T] 
holds true. 

Of course, the backward equation itself is ill-posed, as even if there is exis- 
tence and unicity of the solution (e.g. if the final condition u{T) comes from 
a resolution of the forward equation over the same time period), it does not 
depend in a continuous way of the data. 
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The proof is straightforward by using the following Cole-Hopf transforma- 
tions |Col51l [Hop50l : 



(73) 



Q / U{t,s) ds 

-2v—, v{t,x) ^ v{t,0)e Jo 



V 



Q 7j / U{t,s) ds 

u = -2v-^, v{t,x) =v{t,0)e Jo 

V 

in the forward and backward equations respectively. These transformations al- 
low us to consider the same forward and backward problem, but on the heat 
equation. The Fourier transform gives the existence and uniqueness of a solu- 
tion to the forward and backward heat equation, and the equality between the 
forward v and backward v solutions. Equations ([73]) extend the result to the 
viscous Biirgers' equation. 



4 Non viscous transport equations 

4.1 Linear case: proof of theorem [3] 

In this section we prove theorem [3l 

The first two points of the theorem are easily proven as in theorem [T] with a 
vanishing viscosity. 

Thus we only prove the third point. To do so, we recall that the curves 
(s, tjj{s, x)) are the characteristics of the direct equation ([8]-F) with K = 0, such 
that {s,'4j{s,x))\s=o = (0,x) (see |CH62| lEva98| for characteristics theory). 

For the forward equation ([Sl-F), this change of variable gives 

dsw{s, ip{s, x)) = -if (V'(s, x))w{s, ip{s, x)) (74) 

So that 

w{s,ip{s, x)) = 'w{0, x) exp ^ — I K{'il}{a,x)) d(T^ (75) 
And in particular for 'w{T) we have 







w{T, il>{T, x)) = w(0, x) exp ^- K{^j{a, x)) daj (76) 

For w we have similarly 

d,w{s, V'(s, x)) = K'{i;{s, x))w{s, ij{s, x)) (77) 

So that we have: 

I 

w{s,'ip{s, x)) = w{T,^jj{T,x)) exp \ — / K' {'ip[a,x)) da 



w{T,^{T,x)) exp — / K' {ip{a, x)) da 



(78) 
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Using and ([TS]) we get 

w{s,ij{s,x)) 

= w{0,x) exp — / K{ip{a,x)) da exp — / K'{'ijj{a,x)) da 



= w{s,ijj{s,x)) exp K{^p{a, x)) da^ exp J K {ij,' {a, x)) da^ exp J K'{ilj{a,x))da 
= w{s,^p{s, x)) exp J K{ip{a,x)) + K'(Tp{a,x)) da^ 

(79) 

4.2 Non linear case: proof of theorem [4] and proposition 

m 

From equation (jlSp . we deduce that the forward error if satisfies the following 
equation: 

dfW + wdxW + UobsdxW + wdxUobs = -Kw (80) 
By multiplying by w and integrating over f2, we obtain 

i^t ( / ) + / w'^dxW+ / {uobswdxW + w'^dxUobs) = - / (81) 

Some integrations by part give the following: 

dt{\\w{t)f)^ [ {-2K - dxUobs)w^ (82) 

We set M = \\dxUobs\\ oo, and as K does not depend on x, 

mw{t)\\^) < {-2K + M)\\w{t)\\' (83) 
We have a similar result for the backward error: 

d,i\\mr)<i-'2K'+M)\\wit)r (84) 

We first consider the first point of theorem HI i.e. K{t, x) = K. Gronwall's 
lemma between times t and T gives 

lk(r)f < e(-2^^+^^)(^-*)||u;(t)f (85) 
\\w{tW < e(-2^^'+^^)(^-*)||55(T)f (86) 

from which equation is easily deduced. 

In the second case, i.e. Kit^x) = Kl^t^ i^^[t) and by successively applying 
Gronwall's lemma between times and ii, ti and t2, and t2 and T, one obtains 
equation ([TI 
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Finally, in the case K{t, x) = K{x), by considering a similar approach as in 
sectional, i.e. using the characteristics of the direct equation p^ F) (resp. B), 
it is straightforward to prove that 



from which equation (PO)) is easily deduced. 

5 Conclusion 

Several conclusions can be drawn from all these results. First of all, in many 
situations, the coupled forward-backward problem is well posed, and the nudging 
terms allow the solution to be corrected (towards the observation trajectory) 
everywhere and with an exponential convergence. From a numerical point of 
view, these results have been observed in several geophysical situations, and 
many numerical experiments have confirmed the global convergence of the BFN 
algorithm |AB08j . 

The second remark is that the worst situation, i.e. for which there is no 
solution to the BFN problem, is the viscous Biirgers' equation. But in real 
geophysical applications, there is most of the time no theoretical viscosity in 
the equation, and one should consider the inviscid equation instead, for which 
some convergence results are given. From the numerical point of view, these 
phenomenon are easily confirmed, as well as the exponential decrease of the 
error w. But we also noticed that if the observations are not too sparse, the 
algorithm works well even with a quite large viscosity. 

Finally, these results extend the theory of linear observers in automatics 
|Lue66| : instead of considering an infinite time interval (only one forward equa- 
tion but for T +00), one can consider an infinite number of BFN iterations 
on a finite time interval. This is of great interest in almost all real applications, 
for which it is not possible to consider a very large time period. 
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(87) 



and then. 




(88) 
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